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Abstract 

From the asymptotic expansion of the ground state of the Gross-Pitaevskii equa- 
tion in the Thomas-Fermi limit given by Gallo and Pelinovsky [GP], we infer an 
asymptotic expansion of the kinetic, potential and total energy of the ground state. 
In particular, we give a rigorous proof of the expansion of the kinetic energy calculated 
by Dalfovo, Pitaevskii and Stringari [DPS] in the case where the space dimension is 
■ 3. Moreover, we calculate one more term in this expansion, and we generalize the 

CN ■ result to space dimensions 1 and 2. 

in 
o 

1 Introduction 

After recent experiments with Bose-Einstein condensates, new interest has been stimulated 
in the Gross-Pitaevskii equation with a harmonic potential, taken here in its adimensional 
'. form 

ieu t + e 2 Au + (1- \x\ 2 )u - \u\ 2 u = 0, x G R d , teR + , (1.1) 

where the space dimension d is one, two or three, u(t,x) G C describes the wave function 
of a repulsive Bose gas, and e is a small parameter that corresponds to the Thomas-Fermi 
approximation of a nearly compact atomic cloud [FJ |T] . 

A ground state of the Bose-Einstein condensate is a positive, time-independent solution 
u(t,x) = rj E (x) of the Gross-Pitaevskii equation ( II. ip . Namely, 7] £ : M. d t— y M. satisfies the 
stationary Gross-Pitaevskii equation 

s 2 At] £ {x) + (1 - \x\ 2 )r] e (x) - rf e (x) = 0, x G R d , (1.2) 

Tj e {x) > for all x G M d , and r\ e has a finite energy E E (rj E ), where E £ is given by 

E e (u) = [ (s 2 \Vu\ 2 + (\x\ 2 - l)u 2 + \vA dx. 

./rod V 2 / 
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In dimensions one, two and three, provided e is sufficiently small, existence and uniqueness 
of a radial ground state rj £ is known |IM[ IGPj . It is also well known |IM[ IAAB] that i] £ (x) 
converges to r] (x) as e — > for all x G M?, where i] is the so called Thomas-Fermi 
approximation 

(1 - |x| 2 ) 1/2 for \x\ < 1 



^ (X) = \ for|x|>L (L3) 

Several quantities such as the kinetic energy of the ground state 

EkM = £ 2 I \^Ve\ 2 dx, 

can not be accurately approximated just by replacing the ground state r\ £ by 770, because of 
a logarithmic divergence at the boundary region \x\ = 1. In [DPS], in the three-dimensional 
case, Dalfovo, Pitaevskii and Stringari give the correct behavior of the order parameter in 
the boundary region, which for instance provides the first term in the asymptotic expansion 
of Ek j£ (r) e ) close to e = 0. The purpose of this work is to show how we can obtain rigorously 
asymptotic expansions (theoretically at arbitrarily high order of accuracy) of quantities 
such as the kinetic energy of the ground state Ek, £ (rj £ ), its potential energy 



E p (Ve)= I (\X\ 2 - l) V 2 £ dx } 



or its total energy E e (rj e ). In particular, for d — 3, we give a rigorous proof of the expansion 
of the kinetic energy of the ground state given in [DPS], we calculate one more term in this 
expansion, and we generalize the result to one and two-dimensional cases. The calculation 
of these expansions rely on the expansion of r] e into powers of e in the limit £ — Y 0, which 
was established in [GPJ. So, let us first summarize the main ideas which provide this 
expansion of rj £ in |GPj . Since r\ £ is radially symmetric, we can define a function u e on 
J £ :=(-oo,e- 2 / 3 ]by 



1 — \x 



r ]£ (x)=e 1 / 3 v £ [—^Lj, xER d . (1.4) 

Next, we rewrite equation ( II .2 j) in terms of the new real variable y — (1 — \x\ 2 )/e 2 ^ 3 . It is 
equivalent for r\ £ to solve ( 11. 2ft and for v £ to solve the differential equation 

4(1 - e 2 l 3 y)v"M - 2e 2 / 3 du' £ (y) + yu e (y) - v 3 e (y) = 0, y G J £ . (1.5) 

Let iV ^ be an integer. We look for v £ using the form 

N 

v £ (y) = J2e 2n/3 "n(y) + e 2 ^ 3 R N , £ (y), y G J £ . (1.6) 

n=0 

Expansion ( 11.61) provides a solution of equation ( 11. 5ft if {v n }o^n^N and Rn, e satisfy equa- 
tions (JE}, ([LED and below. 



2 



vq solves the Painleve-II equation 



4^(y) + yv (y)-v 3 (y) = 0, y G R, 

for 1 ^ n ^ N, v n solves 

-Au > ;(y)+W (y)u n (y) = F n (y), y G R, 

where 



(1.7) 



and 



W (y) = 3u 2 (y) - y 



F n{y) = - Yl Vn^y^nziyWAy) - 2di/ n _ 1 (y) - ±yi%_ x (y), 



ni + n 2 + n 3 — n 



Rm F solves 



- 4(1 - e 2/3 y)R" Ne + 2e 2/3 dR' Ne + W R N ,e = F N>e (y, R N>e ), y G J £ 



;i.9) 



where 



2JV-1 



F N ^(y,R) = -(4y^ + 2dv' N )- 



2n/3 



E 



"ni "713 



/ 



2N 



3^2n/3 



v 



2 

/ 



2JV+1 



fl- 3 ^ e 2n / 3 u n . (N+1) \R 2 -e< N+ ^ 3 R 3 . 



n=N+l 



Notice that for ^ n ^ N, u n (y) is defined for all y G R and does not depend on e, 
whereas R^^{y) is only defined for y G J £ . 

In order to describe accurately the convergence of r] £ to 770 by the first term (corre- 
sponding to n = 0) in expansion (ll.Op . z/ shall be chosen in such a way that 



e xj3 v, 



\-x t 



y 1 — |x| 2 if|x| < 1 
if Id > 1, 



which means that uq has to satisfy the asymptotic behaviour 

u o{y) ~ y 1 ^ 2 as V — > +°° an d ^0(2/) ~ ^ as y — >■ —00. (i-io) 

The Painleve-II equation is known to have a unique solution v$ with the asymptotic be- 
haviour (11.1 Oft . This is the so-called Hastings-McLeod solution. Moreover, the behaviour 
of v Q (y) as y — > ±00 has been studied in details, for instance in |HM] . [M] , |FIKN] . Some 
of its properties are summarized in the next proposition. 
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Proposition 1.1 \HM\ [Ml \FIKN^ The Painleve-II equation fli.Tp admits a unique solu- 
tion i> e C°°(R) which satisfies U.10\) . This solution z/ is strictly increasing on R. The 
behaviour of v as y — )■ — oo zs described by 



whereas as y — » +00, 



00 j. 



where bo = 1, b\ = 0, and for n ^ 0, 



g ra+l j n n+1— i 

6 n+2 = 4(9n 2 - l)b n - 2 b 'n b n+2-m - - ^2 b l b m b n+2-l-m- 

m=l i=l m=l 

Once ^0 has been chosen to be the Hastings-McLeod solution of the Painleve-II equation 
(jl.7p . we show by induction on n that there is a unique way to construct the sequence 
(y n )n^i C if°°(R) such that (11. 8p is satisfied for every n ^ 1. Moreover, like for u , the 
asymptotic behaviour of the z/ n 's can be precisely described, as it is shown in the next 
proposition. 

Proposition 1.2 IGP\/ For every n ^ 1, 

00 

v n{y) ~ zA 2 ™ ^ 3n,mf 3m/2 /or some {# n ,m}meN, 

j/— >+oo ' 
m=0 

and i/ n (y) « 0, 

y— >— oo 



where (3 



-5/2 ifd=l, 
1/2 i/d=2,3. 



Remark 1.3 T7ie coefficients g n ,m can be calculated explicitely by plugging this expansion 
of v n into ( fl.gj) . For instance, V\{y) ~ H^~ d ) y-9/2 if d = \ } whereas V\{y) 

l^y-V 2 if d = 2,3. 



J/— S> + 00 "* 3/— >+oo 



Finally, we close the argument by constructing a remainder term -Rtv )£ that solves 
equation (11.91) . We also prove suitable estimates on Rn, £ which ensure that the last term 
in the right hand side of (11. 6ft is indeed small compared to the other ones, and that the 
solution of the stationary Gross-Pitaevskii equation ( II ,2p that we have constructed by 
(11.41) . ( II. 6p and our choices of the ^„'s and R^ >£ is indeed the unique ground state of ( II ,2p 
(in particular, it is positive). More precisely, we have the following result. 



4 



Proposition 1.4 !GP\J For every N ^ 0, there exists > and > such that for 
every < e < En, there is a solution R^ :£ G C°° fl L°°(J £ ) of equation U.9\) with 



\Rn,e\\l™(j £ ) ^ C 



N, 



-2/3 



R N M 2 W,{y){l - e^yy/^dy ^ C N 
1 — \x\ 



e H 



2/TTD d\ 



and S N)E : x i-> R N , £ . ^ 2/3 
such that the unique radially symmetric ground state of equation U.fy) in L 2 (IR d ) writes 

A' 



T] £ (x) = 6 



n=0 



1 — |x| 

^2/3 



^V/3+l^ 



N.e 



1 — |x| 

p2/3 



x G 



1.13) 



Remark 1.5 The estimate on the L°° norm of R N ^ £ written in [GP] was H-Ra^Hl 00 ^) ^ 
Cjv£ _ ^ -1 ^ 3 - The estimate we have in Proposition \1.4 above is a direct consequence of this 
inequality with R Ne replaced by Rn+i,e, taking into account that d ^ 3, Rm, £ — v N+i + 
e 2 ^RN+i, £ o,nd z/tv+i G L°°(M). The other estimate on Rn, £ is a byproduct of the proof of 
Proposition \l-4\ given in jGPf . Indeed, Rn,e is obtained there thanks to a fix point argument 
in a ball B H i(R% e ,Ce 2 ^ 3 ), where \\R% £ \\h^ < 1 and 



u 



Hi 



-2/3 



(4(1 - e^y)^\uf + (1 - E^ y y/^W,{y)u{ y y) dy. 



The asymptotic expansion of r\ £ given by (II .4ft - (11. 6ft as well as the precise description 
of the behaviour of v n {y) as y — > ±oo (for n ^ 0) given in Propositions 11.11 and II. 2| 
enable us to calculate expansions for kinetic, potential and full energy of the ground state. 
Concerning the full energy, we get the following expansion. 



Theorem 1.6 For d = 1, 

2 



Es(Ve 



e 2 \ne + 

15 3 



1 1-tH-i i 

dt - - 

t 2 



My) - y + - x {^i} ) d v 



+ 

For d = 2, 
EM 



1 1 f +OQ ( 2 \ /" H 



MyfMy)dy 



e 8 / 3 + 0(e 3 ). 



7T 2lT 



-2tt 



e Ine 



71 

2 Jo 



My) My) + o 1 ^} dy 



My) - V + -^i} ) dy + n 
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For d = 3, 



+ 



2vr 



+ 



167T 


47T 


105 


T 




-*) 


Jo 


t 


p+oo 





e lne 



-dt-ix J I u (y) 4 - y 2 + ~l{^i} J + ~ 



7T 



n / 2/ ( Mv) -V+ + _1 {v>i} ^ - 4vr / (^o(2/) 3 ^i(2/) + Ifo^o}) ^ 

^ j-oo v y j J -oo 



Remark 1.7 iVoie £/ict£ the notation v\ does not represent the same function in the three 
expansions given in Theorem ] 1.61 Indeed, the equation U.8\) satisfied by v x depends on the 
dimension d through F\. 

The rest of the paper is organized as follows. In section [2] we calculate asymptotic 
expansions of E £ (i] £ ) and prove Theorem 11.61 In section [3] we calculate the asymptotic 
expansion of the potential energy E p {j] e ). In section HJ we deduce the expansion of the 
kinetic energy from the results of the two previous section, and we rediscover the expansion 
found by Dalfovo, Pitaevskii and Stringari in [DP S] on a formal level. In the appendix, we 
prove a key lemma which is used on many occasions in the calculation. 



2 Expansion of E £ (rj £ ) 

We are interested here in the behaviour of E e (r) e ) as e — > 0. First, if we multiply (II. 2p by 
r) £ and sum over W 1 , we get 

Ee(Ve) = -~ / Veixfdx. (2.1) 

From the convergence of r\ £ to t/q in L p (M d ) (for any p e [1, +00]) as e — > [GPj . we already 
know 

1 , ( -8/15 if d= 1 

E £ (Ve) — t -o / Vo(x) A dx = { -tt/6 if d = 2 
ZjRd [ -16tt/105 if d= 3. 

Next, we calculate some correction terms in the asymptotic expansion of E e (i] e ) as 
e ->■ 0. From ([2TT]) . fTOj) and (OJ) we infer 
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Ee(Ve) + 
1 

~2 



1 



where 



r)o(x) 4 dx 

d 

(r] £ (x) 4 - r] (x) 4 ) dx 
1 — \x 



dx 



+ 



r 4/3 



l-r" 

e 2 / 3 



2\ 4 



1 — r ; 



"T 

£! 

" 4 
" 4 



3 d— 1 1 



3 d— 1 1 



-2/3 



-2/3 



My) 4 - V^ 4 ) (i - £ 2/3 y) d/2 -^ 

((^(y) + e^R^MY - v^ 4 ) (1 " e^yY^dy 



h 
h 
h 
h 
h 
h 



2 



= £ 



j —oo y 

f (My) 4 -V^ 4 + 2 -Mv>,i}) (i-^ 2/3 y) 
j —oo \ y / 

— 2/3 

4£ 2/3 T ^ (y) 3 it:o, £ (y)(i- £ 2/3 y) d/2 -Vy, 

J — oo 

— 2/3 

6£ 4/3 f " ^ (y) 2 it:o,e(y) 2 (i - £ 2/3 y) d/2 -V y , 

«/ — oo 

— 2/3 

4£ 2 f ^o(y)it:o,e(y) 3 (i - £ 2/3 y) d/2 -My, 

J —CO 

— 2/3 

f it!o, £ (y) 4 (l - e^y) d '^dy. 



2/3,Ad/2-l 



Next, we give an asymptotic expansion as e — > of each of the I/s, with a 0(£) 
remainder. For this purpose, the following lemma, which is proved in the appendix, will 
be convenient. 
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Lemma 2.1 Let g : R i— >• R fre a bounded function, such that g(y) = 0(exp(y)), and 

y— > — oo 

(7(3/) = 0(y~ a ), where a G R. T/ien 

y— >+oo 

gtoO-ffiy^dv (2.2) 
O^- 1 / 3 ) if a ^1/2 

l^g(y)dy + 0(e 1 / 3 ) ifa^z/2 (2.3) 

»(y)dy - (I - 1) ^ 2/3 J-T w(v)dv + « > 5 /2. 

Expansion of Jx. From the change of variable i = e 2 ^ 3 y, we get 

— 2/3 

Jj = _ 2 f (l_ £ 2/3 y )d/2-l^y 

7i y 

x dt 



-2 / (l-t) d / 2 



= _ 2 r LJ-i'-e'-i^j-n-r'-i, 

J £ 2/3 t J t J t 

= |l» e - 2 jf (1 - t) 7"" 1 ^ - 2 g - l) <*» + O(^). 

Expansion of /2- We apply Lemma [2.11 to the function 

2 

#o(y) := ^o(y) 4 - Vv+ 4 + 
Note that Proposition 11.11 yields #o(?/) = C(exp(y)), and po(y) — 0(y~ 5 / 2 ). Thus, 

J/— )■ — oo j/— )- + oo 

/e" 2 ' 73 r+oo / j \ /.+oo 

g (y)(l - e^yfl^dy = / g (y)dy - - - 1 £ 2 / 3 / j^dy + 0(e) 
-oo J —oo \* / J —oo 

Expansion of J 3 . From fll.6p . we get 

— 2/3 

»/ — oo 

- e -3/3 / k \ 

= ^ \ u (yf $> 2 ^Vi(y) + ^ k+1 ^R k+1 ,M (1 - e^y) d ^dy. 
J -°° \j=o J 

From Propositions 11.11 and 11.21 we have 



My) = y 1/2 -ly- 5/2 + o(y- 4 ), 

y^+oo I 



f V^ 3/2 + 0(z/- 3 ) if d = 2,3 

^^^ooi O^ 9 / 2 ) if d=l, 



and fo(y), vi(y) ~ 0. In particular, for d = 2, 3, 

J/-S-+0O 



n>(v)Vi(i/) = (y 1/2 + 0(r 5/2 )) 3 ( V^^ 3/2 + °^ 3 ) ) = + °^ 3/2 )' 

a result which also holds for d = 1. Thus, Lemma [2.11 implies 

— 2/3 

4e 2/3 f i^ViCyXl-e 2 ^)^!/ 

«/ — oo 

— 2/3 —2/3 

= 4, 2 / 3 T" g 1 {y){l-e 2 "y) d ' 2 - 1 dy- ^ f ' ±=±{1 - e^yf^dy 

J-oo JO Z 

— 2/3 



-oo 

-+oo 



d 



4^ + 4e 2 / 3 / fc(y)<fo + 0(e), (2.4) 



d 

where 

gi(y) = MvTMv) + ^-Mvm = °(^ 3/2 )- 

Z y^+oo 

Next Proposition ll.2l provides. for j ^ 1, i/, +1 (y) = 0(y _7//2 ), and therefore z/ (?/) 3 z/ ?+1 (|/) 

y-^+oo 

0(y~ 2 ). Thus, Lemma [2.11 implies 

— 2/3 

4 £ 2 ( 1 +iV 3 f ^(yJVdlKl-^y)^ = 0(e 2 ( 1+ ^ 3 ) = 0(e 4 / 3 ). (2.5) 

J — oo 

Finally, if £ ^ 1, 

4e 2/3 / i^V^^iWl/Xl - e^y f^dy 

— 2/3 

< 4^ +2 )/ 3 iii? fe+1 , £ ii LO o (j£) r z, (y) 3 (i-£ 2/3 2 /) d/2 - 1 ^ 

«/ — oo 

/ pO re- 2 / 3 ^ 

< £ 2( fe+ 2)/3 / ^ (y)8(1 + lyl) l/2 dy + ^-2/3)3/2 / (1 _ ^ y) ^ dy 

\J -oo JO 

< e 2(*+2)/3(l + e -5/3 ) = 0(e)j ( 26) 

provided k ^ 2. With such a choice of fc, the combination of estimates ( 12. 4p . ( 12. 5 p and 
(1231) yields 



1 — d f + °° 
h = 4— p + 4e 2/3 J gi (y)dy + 0(e 



Estimates on /4,/ 5 ,/ 6 . As mentionned in Proposition Vq is an increasing function 
on R, uo(y) — > and Uo(y) ~ y 1 ^ 2 - On the other side, it is proved in |GPj that 

Wo(y) — 3^(2/) 2 — y > for y G R. Thus, there exists C > such that for every j/ El, 

max(z/ (y) 2 , z/ (y), 1) CW (y). 

Thus, it follows from the estimates on R e stated in Theorem 11.41 that 

J 4 = 0(e 4 / 3 ), h = 0(e 2 ) } h = 0(e^). 

Combining the asymptotic expansions of all of the I/s, we obtain finally 



E e{Ve) + I / rj (xfdx 



4 1 1 



4, ( A i -d n r 1 (1 -ty/ 2 - 1 - 1 , 

-ln£+ 4 — 2/ ^ '- dt + 



d 







t 



9o{y)dy 



j r+oc r+00 \ 

1 - -)(2 + y yy (y)dy) + 4 y ^i(y)dy) £ 2/3 + 0(e) 



3 Expansion of E p (r} £ ) 



1 IU 114 



A calculation similar to the one we made to compute E e (r) e ) + oWVo 



gives 



E p(Ve) ~ I {\ x \ 2 - l) r] (x) 2 dx 

y ((My) + e^RoM) 2 ~ ^) (1 " e^yf^dy 



3 d— 1 1 



-2/3 



"IE* 



where 



— 2/3 

Ji = f --l^Cl-^y)^" 1 ^, 

J — 00 y 



J. 



00 

£ -2/3 



J, = 2 £ 2 / 3 



J 4 = £ 4/3 



y(^o(y) 2 - y+) + ~i {2/> i } J (1 - e 2/3 y) d/2 -My, 
-2/3 

yz/ (y) J R , £ (y)(l-e 2 / 3 y) d / 2 -My, 
yi^(y) 2 (l-e 2 / s y) d / 2 - 1 dy. 



00 
-2/3 



The method used to calculate asymptotic expansions of the J/s for j = 1, 2, 3, 4 is similar 
to the one we used for the I/s. More precisely, 
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Expansion of J x 



1 2 

Ji = -L = - me 

2 3 



•e-'^-'-'a-d-H^+ofr*) 



Expansion of J 2 . From Lemma [2. 1^ 
J 2 = 

where 



/+oo /J \ /.+oo 

fla(y)dl/ - (J - lj £ 2/3 y V92{y)dy + 0(e) 



g 2 (y) = y(My? - v+) + ~Mv>i} = o{ y -^). 

y y->+oo 



Expansion of J3. We infer from Propositions (11. ip and (11.21) that 

1 -d 



yMy) u i(y) 



y^+oo 2 



+ o(r 3/2 ). 



We put 



9s(v) ■= yMy)My) + -^r^ 1 ^} 



2 — ' J/— > + oo 



Q(y 



-3/2N 



Like in the calculation of J 3 , we notice that for every j ^ 1, yvo(y)vj+i(y) 
and we deduce 

1 — d f +OQ 
Js = 2— + 2e 2 ' 3 j g 3 (y)dy + 0(e). 

Expansion of J4. Like in the estimate on J4, we have 

J 4 = 0(e 4 / 3 ). 

As a conclusion, we get 



y— >+oo 



Q(y 



-2\ 



e p(Vs) + / m(%) 4 dx 



jd-l I 



2 n / i-d r 1 (1 - ty/ 2 - 1 - 1 , r +o ° ... 



+ ( (1 - |)(1 + / yg 2 {y)dy) + 2 / g 3 (y)dy ) e 2/3 + 0(e) 
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4 Expansion of E ki£ (r} £ ) 

Let us multiply (11. 2p by r] £ and sum over M. d . We get 

E k>£ {r] £ ) = 2E £ { V£ ) - E p ( V£ ) 

Thus, the asymptotic expansions of E £ {rj e ) and E p {rj £ ) obtained in the previous sections 
ensure 



E k ,e{Ve 



?d-l I 



2, f n l-d f 1 (1 - t Y' 2 ' 1 - 1 , 

-lne+ 2 / - - dt 

3 V d In t 



{%{y) -92(y))dy 



-((1-2)( 1 + 



+oo 



y(9o(y) - 92{y))dy) + / (A 9l (y) - 2g 3 (y))dy e 2 ^ + 0(e 



2 , (c^-d 



1 (1 - tfl 2 - 1 - 1 



o 



+ 



+oo y 

(My) 2 (My) 2 -y) + ~i{ V >i})dy 
oo y 



+ ((!--)(/ (yMy) 2 (My) 2 -v) + W) dy) 



+oo 



((2z/ (y) 2 - y)Mv)My) + (d~ l)l{^o } ) e 2/3 + o( 



+2 



In [DPS], the author study the kinetic energy 

h 2 f 

E kin = — / |W>(r)| 2 (ir 
2m J R 3 

of the ground state if) of the Gross-Pitaevskii equation with an isotropic harmonic trap 



ft 2 a , / \ ( 1 2 i i2\ , / \ ^h 2 a 
— AiP(v) +U- -mu HO \v\ if>(v) - — 



|^(r)| 2 ^(r) = 0, 



(4.1) 



where /i > is a chemical potential, a > is the scattering length. Like in [DPS], we 
denote by R the radius of the condensate, defined by 



2 d2 

H = -mu HO R , 
and we introduce the harmonic oscillator length 

( h y/ 2 

a HO = 

\muHoJ 
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as well as the maximal value a of the wave function ipTF of the Thomas- Fermi approxima- 
tion: 



R 



a 



(8?< a) 1/2 ' 
Then, the total number of particles is 



^ TF (r) = a ( 1 - ^ if |rK 1. 



N 



ip TF (r) 2 dr 



R 5 



r <1 



1 ^)QjQj j-^Q 



The change of variables 

ip(r) = au(r/R), 

maps the ground state ^ of (14. ip to the solution u = r] e of (II. 2p . where 

a HO 



R 2 ' 



Let us now use our expansion of E k>e (j] £ ) to calculate the expansion of E kin . 

h 2 f 

2m J R3 



h 2 R 3 
2m SirajjQa 

h 2 R 3 



ima HO a 



s- 2 E k:£ (r] £ ) 
3 R 2 3 



-dt 



+ 



My) - 92{y))dy + 



5NH 2 
~2mR 2 



duo 1 3 , 3 ' 

In — + In 2 + - 

R 2 2 4 



Myf-yMyf 



Moreover, taking into account the properties of the solution v of the Painleve mentioned 
in Proposition ll.l[ we infer 

i C k - ^ +h f-y-\- C - li r) *■ 

In order to compare our result with the one in (DPS], we introduce the function = 
2~ 1 / 3 z/ (— 2~ 2 / 3 £), which is the solution to 

y?" - ^ - ^ = 0, ¥>(£) t ~ f(0 c ^ 0. 
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Thus, 



5Nh 2 
~2mR 2 

5Nh 2 
~ 2mR 2 
5Nh 2 
2mR 2 



ln^ + 2-^ln2-3 [ + °° 



In — -2 + -ln2 

o>ho 2 



W + l« S _p)^ + o((^) 4 ' s 



+3 lim 



2 r 



2mi? 2 
2mR 2 



In 



X;^4^- 2/3 -z^) + °((^) 4/3 



in 



-2 + h l2 + 3(c + iln2) + o((^) 4/3 ) 



where like in [DPS], we have denoted 

C= lim ( ! + <p'(£) 2 d£- -hi A] - -In 2. 



5 Appendix: Proof of Lemma 12.1 



If a ^ 1/2, we split the integral in three pieces. On the one side, since d ^ 3 and since the 
maps g and y h-> \y\ l ^ 2 g{y) are in L 1 (M_), we have 



\g(y)\ max(2,2e 2 / 3 |y|) 1 /^ = (1). 



Then, there exists C > such that g>(y) ^ C(l + |y|) x / 2 for every y ^ 0, and since d ^ 1, 

/■e- 2 / 3 /2 /- £ - 2 / 3 /2 

/ ^(l-^y^dy < / C7(l + |y|)- 1 /»2 1 / a dy = 0(e- 1 /3). 
Jo Jo 



Finally, 



-2/3 



£-2/3/2 



^( t / £ 2/3 )(1 _ t) d/2- 



1/2 



£ 2/3 



< Ce 1 / 3 f t - 1 / 2 (l-t) d / 2 - 1 ^ = 0(s- 1 / 3 ). 

A/2 
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If a ^ 5/2, using the Taylor formula, 

— 2/3 

f g{y){l-e 2/ *y) d/2 - l dy 

J — oo 

/+oo i>e- 2 /' A /2 
9(y)dy+ / g{y) ({1 - e 2 "y) d ' 2 - 1 - l) dy 
-oo J —oo 

re~ 2 / s r+oo 

+ / g^a-e^yY'^dy- g{y)dy 

j£-2/3/2 j£-2/3/2 

/+oo [>e- 2 l' A /2 
g(y)dy+ / 
-oo J —oo 

where £ ejJ/ G [0,e 2/,3 y) or £ £jy G (e 2 / 3 y,0), depending on the sign of y. In particular, if 
y s= £" 2/3 /2, & tV < 1/2 and since d/2 - 3 < 0, (1 - 6,y) d/2 " 3 < 2 3 ~ d/2 - Additionnally, 



"+00 

/e-2/3/2 

-(^ - i^ /3 y + - 1)(^ - 2)i(i - ^"^V ) dy 



f 12 g(y)y 2 dy ^ f g(y)y 2 dy + f " C(l + y)- 5 / 2 y 2 dj, = O^ 1 / 3 ) 

J —oo J —oo Jo 

Moreover, there exists a positive constant C such that for t G [1/2,1], g(t/e 2 ^) ^ 
(e 2 / 3 /t) 5 / 2 . Thus, 



/ 3 /2 



-2/3 



A 

g(y)dy-q-l)e^ 



"2/3/2 



OO "J —oo 

+oo 7 /"+00 

- - 1) £ 2 / 3 

2 ; 



g{y)dy 



oo 



yg(y)dy + 0{s) 



yg{y)dy + 0(e). 



oo 



If a ^ 3/2, we start the calculation like in the case a ^ 5/2. Then, again thanks to the 
Taylor formula, 



/ g(y) ((i - s 2 ^ 2 - 1 - 1) dy 

J —oo 



2/3 / 2 



22-^/2^.2/3 



f!-i 

2 



£-2/3/2 



|<7(y)IMdy = 0( £ 1/3 ), 



for some G [0,e 2 / 3 y) or £ £)2/ G (e 2/3 ?/,0). Then, 

f - s 2/3 y) d/2 -Vy = f g{t/e^){l - tf 2 - 1 ^ = 0(e^). 

A-2/3/2 A/2 £ 1 



Finally, 



+oo 



'e-2/3/2 

which gives the result in the case a ^ 3/2. 



g(y)dy = 0(e 1 "), 
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